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(54) Title: RETAIL SHELF INVENTORY SYSTEM 



(57) Abstract 

A retail shelf inventory management system includes memory 
for storing inventory data. Th e inventory data includes a study period 
and a retaiT opening time and a retail closing timefo r eachj^ay wittan 



[SEES^_Spd_for 

each of a number of different inventory items or products. The histori- 
,car&ventory "data4ndp a cust omer service le- 

vel, a proportion of business done for each day within the study peri- 
od, and a variability of demand. An optimizing program and process- 
ing unit are operativelv_ connected J to Ae memory devices for optimiz- 
ing ashelf capatityJ204) for anv-o neach of saidj uIrSluXiLu^ 



itemsTBoth"user selected shelf capacity values and the optimized shelf 
"capacity value can be separately utilized for calculating target stock le- 
vels (205), order-to-levels, estimated sales estimated lost sales (208), 
average inventory values (210) and customer service levels achieved. 
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to provide an improved retail shelf inventory management 
system capable of identifying an optimum shelf capacity on 
an individual product basis; to provide such a system capa- 
ble o^f-rd^ffETfying target stock level s and order-to-levels 
5 separately utilizing the identified optimum shelf capacity 
and/or a user supplied shelf capacity value; and to provide 
such a system capable of identify ingles timated sales and 
customer service levels achieved. 

^ — . 7 

In brief, the objects and advantages of the pre- 

10 sent invention are achieved by a retail shelf inventory 
management system including memory for storing i nventory 
data. The inventory data includes a study peri od and a re- 
tail opening time and a retail closing time foreacji day 
within the^gt udy p eriod anSThistorical inventory data for 

15 the study period for each of a number of different inven- 
tory items or products. The historical inventory data in- 
cluded a movement value, a customer service level, a pro- 
portion of business done for each day within the study pe- 
riod, and a va riability of de mand. Optimizing means are 

20 opera tively connected to said mem ory means for o ptimizing a 
shelf capacit^ for a ny or each of said plurality of in- 
ventory items,.-. Both use r: selected s helf capacity values 
and the optimized shelf capacity value can be separately 
utilized for calculating target stock levels, order-to- 

25 levels, estimated sales, estimated lost sales, average 
inventory values and customer service levels achieved. 

Brief Description pf the Drawings 

The present invention together* with the above and 
other objects and advantages may best be understood from 
30 the following detailed description of the embodiment of the 
invention illustrated in the drawings, wherein: 

FIG. 1 is a block diagram representation of a re- 
tail shelf inventory system according to the present inven- 
tion; and 



WO 90/09638 _ 3 _ PCF/18S9O/O0652 

FIGS. 2-7 are flow charts illustrating the logi- 
cal steps performed by the retail shelf inventory system of 
FIG. 1. 

Detailed Description of the Preferred Rrnhnrfin^ 
5 Referring now to the drawings. In FIG.vl there is 

illustrated a block diagram representation ots retail 
shelf inventory system generally designated by t&e refer- 
ence character 10. 

Appendix I below sets forth a detailed descrip- 

10 tion of a stochastic inventory model according tq the pre- 
sent invention. Appendix I provides a problem definition 
and provides the results by the stochastic inventory |iiodel 
of the invention and includes both a mathematical descrip- 
tion and a corresponding exemplary program specification 

15 written in Pascal for implementing the mathematical de- 
scription, il , ; • 

As illustrated in FIG. l, the retail shelf inven- 
tory system 10 includes a central processing unit 12 and an 
associated memory generally designated by the. reference 

20 character 14. As shown , the memory 14 includes program 

memory space 16 for storing the retail ^helfi inventory pro- 
gram of the invention and a plurality c# data memory Spaces 
18 for storing multiple data files (l)?- i {Q) including both 
user supplied data and computed results data of $ata sets 

25 and arrays. A keyboard 20 is coupled to the central pro- 
cessing unit 12 for entering user selections an* data. A 
display 22 and a printer 24 are also cpupled to the central 
processing unit 12 for reporting .calculated results to the 
user. The retail shelf inventory system 10 can be impler 

30 mented by a personal computer system, for example, such as, 
an IBM PS/2 Model 80 with an associated disk drive 70J4 
bytes memory, or various other commercially available mi- 
crocomputer based systems. 

As set forth below in Appendix I undex^theadjing 

35 M B. Results from the model", the retail shelf inventory 

system 10 provides an optimum shelf capacity on an individ- 
ual product basis. Optimum stocking times when a shelf 
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should be replenished also can be provided by the retail 
shelf inventory system 10. For both the calculated optimum 
shelf capacity and each of the user entered shelf capacity 
values, the inventory system 10 generated results can in- 
5 elude order-to-quantities for reordering, estimated lost 

sales and average inventory stock values for both the store 
and shelf stock as seen by a customer. In accordance with 
a feature of the invention, the daily opening and closing 
store time s and historical values of the proportlons~o f 

10 business done by day within a study period are supplied by 
the user and utilized for generating the user selected re- 
sults to be reported for a particular product. The use of 
actual opening and closing store times together with the 
proportional business done values simplifies mathematical 

15 calculations and can provide increased accuracy for gener- 
ated results. 

Referring now to FIG. 2, there is shown a flow 
chart generally illustrating the overall program logic se- 
quence in accordance with the principles of the invention 

20 for the shelf inventory system 10. The program starts with 
an INITIALIZATION routine performed indicated at a block 
200. The INITIALIZATION routine is illustrated and 
described with respect to FIG. 3 and in Appendix I under 
the heading "Solution: Mathematical Description", 

25 subheadings "l. Notation" and "2. Preliminaries"; under 

the heading "Solution: Program Specifications," subheadings 
"1. Variables" and "2. Initialisation." 

Then an optional STOCKING TIME routine indicated 
at a block 202 indicated at a block 204 shown in dotted 

30 line can be performed to provide optimum stocking times 
when a shelf should be replenished. At least an initial 
stocking time is supplied by the user, then subsequent 
optimal stocking times are calculated to space the subse- 
quent stockings evenly in shop time over the study period. 

35 Subsequent optimal stock times are identified by calculat- 
ing an inverse function of clock as defined in Appendix I 
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under the heading "Solution: Mathematical Description," 
subheading "3. Optimal stocking times;" 

Next an optimum SHELF CAPACITY routine indicated 
at a block 204 is performed for providing an optimum shelf 
5 capacity value. The optimum SHELF CAPACITY routine is 
illustrated and described with respect to FIG. 4 and in 
Appendix I under the heading "Solution? Mathematical 
Description," subheading "4. Calculating optimal. -shelf 
capacity" and under the heading "Solution: Program 
10 Specifications," subheading "3. Calculating shelf 
capacity." 

Next a TARGET STOCK routine indicated at a block 
206 is performed for calculating target stock levels and 
order-to-levels . The TARGET STOCK routine is illustrated 

15 and described with respect to FIG. 5 and in Appendix I 
under the heading "Solution: Mathematical Description, " 
subheading "5. Calculating target stock levels,.* and order- 
to levels" and under the heading "Solution: Program , 
Specifications," subheading "4. Calculating target stock 

20 levels, and order-to levels." 

Next a LOST SALES routine indicated at -a block 
208 is performed for calculating the expected lost sales 
and the customer service level achieved for the study 
period. The LOST SALES routine is illustrated and 

25 described with respect to FIG. 6 and in Appendix I under 
the heading "Solution: Mathematical Description/" 
subheading "6. Lost sales" and tinder the heading ''Solution: 
Program Specifications," subheading "5. Lost sales." 

Next an AVERAGE INVENTORY routine indicated at a 

30 block 210 is performed for calculating the expected average 
shelf inventory seen by the customer, the expected 1 average 
store inventory at the close of business on a particular 
day within the study period and the expected total store 
amount. The AVERAGE INVENTORY routine is illustrated ,and 

35 described with respect to FIG. 7 and in Appendix. I uprier 

the heading "Solution: Mathematical Description subhead- 
ing "7. Average inventory calculations" and und&r the 
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heading "Solution: Program Specifications," subheading "6. 
Average inventory calculations." 

Referring to FIG. 3, there is shown a flow chart 
illustrating the sequential steps performed during the IMI*- 
5 TILIZATION routine. The INITILIZATION routine begins with 
the definition of initial and default values for constants, 
variables, and functions utili2ing user supplied data 
structure definitions. The user supplied data is set forth 
in Appendix I under the heading "A. factors influencing the 

10 model." The user supplied data is included within the list 
of variables used in the program in Appendix I under the 
heading "Solution: Program Specifications, subheading "1. 
Variables" as defined under the heading "Solution : 
Mathematical Description", subheadings "1. Notation" and 

15 "2, Preliminaries." The user supplied data is read and 

copied into the defined data structures, then precalcula- 
tions are performed on the user supplied data on an indi- 
vidual product basis indicated at a block 300. 

As illustrated in Appendix I under the heading 

20 "Solution: Program Specifications, subheading "2. Initial- 
isation," the user supplied time schedule values are con- 
verted from real time into shop time and the proportion of 
total business done and the modified mean demand in each 
interstocking period is calculated indicated at a block 302 

25 and described with reference to FIG. 3 as follows. 

First an order shop time array and a delivery 
shop time array is calculated utilizing the user supplied 
store opening and-closi nq time for each da vJ=0 , . . . , L-l 
IiTTEe^study period and user supplied order and delivery ^ 

<— — - . — — — ■ 

30 real time schedule indicated at a block 304. Next when 
stocking times are user supplied, the user supplied 
stocking times are converted to shop time to generate a 
stocking shop time array indicated at a block 306. Other- 
wise when only an initial stocking time is user supplied, 

35 then the stocking shop time array is calculated indicated 

at the block 306. Next the proportion of business done and 
the modified mean demand during each stocking interval is 
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calculated indicated at a block 308.. Error checking 
routines indicated at blocks 310, 312 are sequentially 
performed on the calculated results for each df the above 
described steps. When errors are not detected at the block 
5 310, the results are stored indicated at a block 314 in 
predefined corresponding data structure for results of 
supply, stock and demand data sets. 

Referring now to Fig. 4, there is shown a flow 
chart illustrating the logical steps of the SHELF CAPACITY 

10 routine performed for calculating an optimal siielf capac- 
ity. Appendix I under the heading "Solution: Mathematical 
Description", subheading "4. Calculating optimal shelf ca- 
pacity" and under the heading "Solution: Program Specifica- 
tions, subheading "3. Calculating shelf capacity" describes 

15 the SHELF CAPACITY routine in detail. The sequential steps 
begin by identifying one of Methods 1 indicated at a block 
400, Method 2 indicated at a block 402 or Method 3 
indicated at a block 404 to be utilized for calculating the 
optimal shelf capacity. A particular Method 1, 2, or 3 is 

20 performed responsive to a user entered method . selection for 
the optimum shelf capacity calculation. Next for either of 
the Methods 1, 2, and 3, historical demand inventory data 
for the particular product, the proportion of business 
done, the modified mean demand during .each interstocking 

25 interval and the desired customer service level is col- 
lected indicated respectively at block 406, 408 or 41,0. 

When Method 1 is selected a root finding routine 
is performed to identify the optimum shelf capacity C. 
This root finding routine evaluates repeatedly the expected 

30 overall lost sales indicated at a block 412 for the entire 
study period utilizing user supplied data,, ; the calculated 
values of proportions of business done and modified mean 
demand for each interstocking period, the demand and the 
customer service level to be achieved as set forth in Ap- 

35 pendix I under the subheadings "Method 1»" If the root 

finding fails indicated at a block 414, then an e^jror flag 
is set and no value is stored. Otfierwise the result from 



WO 90/09638 _ PCT/US90/00652 



the root finding calculation is stored indicated at a block 
416 as the optimum shelf capacity. 

When Method 2 is utilized for finding the optimum 
shelf capacity C, then an index k is found for which the 
5 proportion of business done in the j th interstocking period 
is greatest. Then the optimum shelf capacity is identified 
to achieve the desired customer service level during the 
identified busiest interstocking period indicated at a 
block 418 rather than the entire study period utilized in 

10 Method 1. The Method 2 routine is set forth in Appendix I 
under the subheadings "Method 2." When an error indicated 
at a block 420 does not result from the calculation of the 
optimum shelf capacity of Method 2, then the computed 
optimum shelf value is stored indicated at a block 422. 

15 Method' 3 identifies an optimum value for shelf 

capacity utilizing a block of the busiest or worst r con- 
secutive days identified in the entire study period 
indicated at a block 424 as set forth in Appendix I under 
the subheadings "Method 3." The worst r consecutive days 

20 can occur within a single or multiple interstocking 

periods. When an error indicated at a block 426 does not 
result from the calculation of the optimum shelf capacity 
of Method 3, then the computed optimum shelf value is 
stored indicated at a block 428. 

25 When the user has supplied one or more values for 

the shelf capacity indicated at a block 430, for example , 
such as, a current shelf capacity and a selected shelf 
capacity then the user defined C value or values are stored 
indicated at a block 432 for subsequent computations. 

30 Referring to FIG. 5, there is shown a floy } chart 

illustrating the TARGET STOCK LEVEL routine. The sequen- 
tial steps begin with identifying first a stored shelf 
value C indicated at a block 500. The identified stored 
shelf value C can be either the calculated optimum shelf 

35 value or one supplied by the user. Next the expected sales 
in the j tt interstocking period are calculated indicated at 
a block 502 using the formula set forth in Appendix I under 



WO 90/09638 -9- PCT/IJS90/00652 

the heading "Solution: A Mathematical Description , 
subheading "Method 3." Then an effective delivery" time is 
calculated indicated at a block 504 as set f or£h in 
paragraph (5i) in Appendix I under the heading "^plution: 
5 Mathematical Description", subheading "5. Calculating 
target stock levels, and order-to levels." Next an 
internal target stock level value W k is calculated 
indicated at a block 506 using the formula set £ pjrth under 
paragraph (5). The calculated internal stock .level is then 

10 stored for subsequent computations. Next an ordsr-to-level 
is calculated indicated at blocks 508, 510 utilising the 
shelf value C and the stored internal target stock level by 
calculating the expected sales between the order time and 
the delivery time using the formulas set forth under para- 

15 graphs (5a), (5b). The order-to-level is identif^d to en- 
sure that the stored internal target stock level value W k 
is satisfied at the delivery time. 

Referring now to FIG. 6, there is shown , a flow 
chart illustrating the logical steps performed for the LOST 

20 SALES routine. The LOST SALES routine is described in Ap- 
pendix I under the heading "Solution: Mathematical Descrip- 
tion," subheading "6. Lost sales" and under the heading 
"Solution: Program Specifications," subheading "5,* Lost 
sales The sequential steps begin with a calculation 

25 indicated at a block 600 of the effective stock level that 
is recursively calculated starting from a selected index k 
for which internal calculated target stock level, is maximal 
utilizing the formulas set forth in paragraphs (§}, (6a), 
and (6b). The calculated effective stock levels are then 

30 stored for subsequent computations. Then the expectejl lost 
sales for the whole study period is calculated indicated at 
a block 602 in accordance with the formula set forth in 
paragraph 6(c). Then the customer service level CSL 
achieved is calculated indicated at a block 604 using the 

35 resultant expected lost sales value utilizing the* formula 
set forth in paragraph 6(d). 
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Referring now to FIG. 7, there is shown a flow 
chart illustrating the sequential steps of the AVERAGE IN- 
VENTORY routine. Appendix I under the heading "Solution: 
Mathematical Description", subheading "7. Average inventory 
5 calculations" and under the heading "Solution: Program 

Specifications, subheading "6. Average inventory calcula- 
tions" describes. the AVERAGE INVENTORY routine in detail. 
The sequential steps begin with calculations indicated at a 
block 700 of the expected amount of backroom stock b and 

10 shelf stock a just after the j th stocking time, utilizing a 
computed expected sales value utilizing the formulas of 
paragraphs 7(a) and 7(b). Next the functions a(t) and b(t) 
are interpolated indicated at a block 702 between the 
calculated a, b values to identify average stock amount on 

15 the shelf and in the backroom in accordance with the 
formulas of paragraphs 7(c) and 7(d). Next numerical 
integration of the a(t) time value is calculated indicated 
at a block 704 as set forth under paragraph (i) "Average 
Amounts Seen by the Customer." Next the average store 

20 inventory is calculated indicated at a block 706 utilizing 
store time as set forth under paragraph (ii). Then a 
numerical integration routine can be performed indicated at 
a block 708 to provide a real time average store inventory 
as set forth under paragraph (iii). 

25 Then the sequential functions or selected func- 

tions of the TARGET STOCK LEVEL, LOST SALES and AVERAGE IN- 
VENTORY routines are repeated for the next stored shelf ca- 
pacity value C. 
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APPENDIX I 
STOCHASTIC INVENTORY MODEL 

Problem Definition si . w 

To derive, on an individual product basis, a stochastic inventory model using the fac- 
tors, given below, influencing sales and hence the stock level and re-otder st^atefgics. 

A, Factors influencing the model 

1. The study period; variable from 1 day to a year. 

2. The daily opening times; for each day in the study period, the opening and closing 
rimes of the store are stated. 

3. The demand cycle; this is derived from the following historical data: 
(i) movement figure over the period 

00 customer service level over the period 

(iii) distribution of business by day (Le» sales percentages by day) 

(iv) variability of demand over the period. 

4. The number of times the shelf is replenished during the study period. 

5. The stocking schedule; each time and day within the study period when the shelf is 
replenished is stated. 

6. The number of deliveries to the store during the study period. 

7. The order and delivery schedule; each time and day when an order is raised, i.e. the 
time at which the amount to be re-ordered is calculated, together with its associated 
time and day of delivery within the study period, 

8. The target stock levels; the total amount of stock in the store after eden delivery. 

9. The actual shelf capacity, the amount of the product currently displayed on the shelf. 

10. The desired customer service level. 
[User definition of factors 
Obligatory 

AltoA4 

First stocking time in A5 
A6toA7,A10. 

Optional 

Rest of stocking schedule in A5 
A8 and A9. 

(If A5 and A8 are not fully defined by the user, then the optimal values axe produced and 
used by the model.) ] 
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B. Results from the model 

1. Optimal stocking times, the times and days during the study period when the shelf 
should be replenished. 

lOnly used if A5 not fully defined by the user.] 

2. Optimal shelf capacity; the amount of the product which should be displayed on the 
shelf to meet the factors A 1 to AS* 

3. For each of the shelf capacities (Le, both Actual and Optimal) the following results 
may be derived: 

(i) Re-ordcr schedule; for each order point a stock level which the store must re- 
order to (hence any stock in the store and any stock currently on order must be 
taken into account when calculating the amount to re-order). 

(ii) Lost sales; the sales not achieved by accepting this shelf size and accounting for 
the backroom. 

(iii) Average inventory; the average amount of stock in the store. 

(iv) Average amount on shelf; the average amount of stock on the shelf as seen by 
the customer. 

4. Actual customer service level; the customer service level actually achieved when 
using the actual shelf capacity (A9) rather than the optimal shelf capacity (B2). 
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Solution: Mathematical Description 



Notation 




L 


length of study period (in days) 


Si 


time of opening on day J>J -0, .... L-\ 


J 


time of closing on day J, J « 0, .... L -1 


r 


historical movement figure for the study period 


ioop% 


desired customer service level 


I00a% 


historical customer service level 




proportion of business done oh day j 


100v% 


variability of demand 


N 


number of stockings of shelf 


T k 


time of stocking of shelf,* =0, .... AM 


M 


number of deliveries 


*>', 


time of)* delivery, y = 0, M A 


Oj 


time of y * order, j=0, M-i 




/• target stock level, 7 = 0, .... AM 


c* 


actual shelf capacity ' J> 



2. Preliminaries 

(i) We input the historical movement figure \l' and the historical customer service 
level a, but what really matters is the demand, m which is the number of times someone 
wants to buy one of the product during the study period. Of course, sometimes a custo- 
mer who wants the product will find the shelf empty and be unable to get the item, but this 
will only happen in a proportion (1 - a) of cases. Otherwise, the customer will take the 
item from the shelf, and the movement figure goes up by one. Hence 

\i' s movement « demand x historical customer service Ifcvel 
so demand n = Ji'/a 

(ii) If the variability of demand is v, the standard deviation wilf be jlv, and so the 
variance of demand, o 2 , will be ftiv) 2 . 

(iii) To calculate the proportions of business done in each intemocking period, we 
make use of the function clock; [Q.L) -> [0,11, which is defined as 

ciock(r) ■ proportion of total business done by time t 
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y=o I v 5 t J 

Wc subsequently extend clock periodically to the whole of the real line, so that 

docket ♦/)■«+ clock(0 

for integers n, and r in [0,L]. The purpose of clock is to convert the real time t into the 
shop time clock(t). As far as the shop is concerned, it is shop time which matters; busi- 
ness proceeds at a constant rate in shop time. 

(iv) We also assume that the study period repeats itself, so that if (7*o,Ti] is the 
zero* interstocking period, the (rV-1)* is (Tn-i> l + 7*al' for example. The proportion of 
business done in the /* interstocking period is simply 

pj m ciock(7; + i) - dockOy) . 

where Ty ■ L +T<j. 

(v) We model the demand in the interstocking period by the random variable 
Yj, which is the positive pan of a normal random variable X), whose variance is pyO 2 and 
whose mean is chosen so that the expected value of Yj, EYj, is equal to the mean 
demand in the period, \ipj. More explicitly, if we define the function A by 

(1) A{z) * 5(i> . 

where 5(z) * P (Z > z) for Z a Af <D. 1) random variable, then 

(2) A(z) = £(Z-z) + ; 
thus we choose \ij to solve 




[Note: The choice of a normal distribution for Z is essentially arbitrary, and any zero* 
mean unit-variance distribution could serve instead. The expression (2) for A stays the 
same, but the explicit formula (I) will change.] 

3. Optimal stocking times 

Suppose we have been told the time To of the zero* stocking, which occurs at the 
shop time to «clock(To). Then the best way to choose subsequent stocking times is to 
space them out evenly through (he study period - but 'evenly' in terms of shop time. Thus 
we define 
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x } » T 0 + J/N , j - 1 AM , 

(with addition mod I) and set the recommended stocking times T\ , ...Jn-\ to be 
Tj m clockinv(t;) , 

where clockinv is the inverse function to clock; clockinv(/) = inf{u: clock (u) > /). This 
function can be expressed explicitly in terms of the 8;, Sj, tj and L very much as clock 
can. Less satisfactorily, it can be evaluated by • general root-finding routine 1 fr&i§€lock. 
The user can, of course, ignore the recommended 7/ and supply his own; , ~- • 

4. Calculating optimal sheir capacity 

If the shelf capacity is C, then assuming that the shelf is full at the beginmng of the 
j * interlocking period, the expected lost sales in the 7* interstocking perioltis given by 
the formula 



(4) X/C) ■ <r1pjA 



which is E{Yj-C)*. / ' ' 

There are now three possible ways of choosing the optimal shelf capacity, depending 
on the user's performance criteria. We consider the first to be the most natural. As usual, 
the user is free to ignore the recommended shelf capacity, and input his own value C* 
instead. 

Method 1. This method should be used if the aim is to ensure that on average 100(J% ot 
the total demand throughout the study period will be met Assuming jhat there is always 
enough available at the start of an interstocking period to fill thd'Sheif^the ejected lost 
sales in the study period is 

/(C) » Z tyC) - 

when the shelf capacity is C, and we simply pick C so that 



/(C) - • 
This will need a root-finding routine. 
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Method 2. This method will ensure that during each inierstocking period, at least 
100(J% of the demand will be met Find the index k for which pj is greatest, and then 
choose C to satisfy 

X*(C) = 0-p)W>* • 

This choice of C will always be higher than in Method 1 . 

Method 3. This method should be used if the aim is to have the shelf large enough to 
carry through the worst r consecutive days in the study period. To do this, let p be the 
maximum of 6/ + • • * + _i as j varies, calculate v as in (3): 



\ip = <Wp A 



f \ 
-v 



and now choose C as in Method 2; 



<rip A 



At this stage, the program has a value of C to work with, either computed by one of 
the three methods, or else input by the user. Now one can compute, for example, the 
expected sales in the /* interlocking period: 

Sj = 5y(C) - £0>aC> = vpj - V,(C) . 

S. Calculating target stock levels, and order-to levels 

The user inputs the number M of deliveries in the study period, the times 

Dot »-J>M-\ at which me y 2X6 made » and lhe timcs °o« — at wnicn * e correspond- 
ing orders arc placed. Some retailers may require to be protected against delays in 
delivery of up to 8, in which case we replace D) by Dj + 8, and since we cannot make use 
of delivered goods until the next stocking after delivery, we immediately replace D) by 
the effective delivery time 

<5i) Dj « inf [T k : T t >Dj) . 



Now if D k = 7„ Dt.t =T, +(l+ 'i (/i 20), our target stock level W t at the delivery 
time D k will be set so as to meet on average the sales to be made in (r,,r, +n+ i]: 
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(5) 



r \ 
Z £(30 a C) 



v C 



Wc take the maximum with C so as to ensure that we start with a full sh§lf at D k * The 
objective is to ensure that just after the Jfc* delivery has been triads the tbwamount in the 
shop (shelf + backroom) is W k (or more). Thus we choose the order-iv level V k to 
achieve this, as follows. 

Firstly we calculate (an upper bound for) the expected sales in (O kl D k \ f ijy summing 
the expected sales 5. for each interlocking interval j contained entirely in (0*,D*L *nd 
then adding a correction for the sales between O* and the next stocking tupe. Explicitly, 
if Ojk occurs in interstocking period m, shop time t after the beginning of the interstocking 
period, then the expected sales between 0* and the end of interstocking period m is 



5(a) S„ - £ 



A C 



Pm 



this is the correction we add in, to arrive at the expected sales U k during (Q k ,D k ]. We 
now define the** order-to level by 

5(b) V* * w k ♦ V* • ' 

The interpretation is this. At time 0*, we have amount jc in the shop (shelf + backroom) 
and we have placed orders for a total amount y which has not yet been delivered; we 
therefore place an order for amount (V* -x -y) + . 



6. Lost sales 

Wc have chosen a shelf capacity C, and target stock levels W 0 , ... Wm-i- Suppose 
that D k * T fl Ojk^i ■ T,+ n +i. Then the lost sales in {D k <D k ¥ \] is exactly 



r «■ n 



The expected value of this is too cumbersome to evaluate explicitly* so we approximate 
the expected lost sales in (D k ,D k ♦ 1 1 by 



(6) A k (C,W k ) = nT/), - 



Wc could now just sum the A*(C, W 4 ) to yield expected lost sales, but this would be 
erroneous for the following reason. Suppose that W\, were all very small, but 
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W 0 was enormous. Then there would be lots of unsold stock at time D\,so effectively 
there would be more in the shop at time D i than W u What we need to do is replace the 
target stock levels W k by the effective target stock levels W k , and then proceed. 

We calculate the effective target stock levels W k as follows. Pick the index k for 
which W k is maximal; suppose for convenience that it is k * 0. Then we set 

6(a) w 0 * Wq , 

and define recursively 

6(b) Wi.t - W t + t v [Wi - M(clock(O i+1 )-Clock(D,)) + A.CC.fe} . 

This is because if we had Wr in the shop at D it the demand in (Pi,D^i) would be 
H{clock(D J+l ) -clock(Dj)), of which on average A.CCW,) would be lost. The amount 
remaining just before D M would be V^-^clcek(D i+ j)-clockCDi)} + A<(C,WU and 
this might exceed VV i+1 . 

Finally, then, we calculate the expected lost sales for the whole study period: 

u-\ 

6(c) Expected lost sales = J ^(C,^) . 

kmO 

From this follows the customer service level achieved: 
6(d) 100 [1 - (Expected lost salesym% . 

7. Average Inventory calculations 

The first step here is to calculate the expected amount on the shelf and in the back- 
room just before and just after the /'* stocking time, Tj. 

Suppose that just after Tj, there is a } on the shelf and bj in the backroom. Just before 
Tj+u there is b, in the backroom and (approximately) -S y )* ■ (ay -£(!>• a C)) + on 
the shelf. Thus 

7(a) bj*x ■ {bj - (^Ad ; )) + 

7(b) a h x - Uflf-Sj)* + bj) a C 

except when 7) is a delivery time D kt in whicvh case bj * aj is simply the effective target 
stock level W k . and aj is just W k a C. 

The second stage is to interpolate between these values to obtain functions a (•), b (•) 
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for the average amount on the shelf and In the backroom. The backroom stock level is 
easy: 

7(c) bit) « bj (t; S t < x h i) . 

For the amount on the shelf, we set ■ 
7(d) fl (f) = (a y -m) v (aj-Sj)* (t y i t < x y >i) . 

[Recall that x } » clock(T,).] From this, we can obtain the average amount on the shelf (or 
in the backroom), however one may choose to define it We suggest some possibilities 
here. 

(i) Average amount seen by the customer. This is simply 

[ l a(t)dt, 

which can be evaluated by a (trapezium rule) numerical integration. 

(ii) Average amount in the store at the end of the day's trading. This amount in the 
shop at the close of trading on day j is simply 

Pj m aCclockty)) + Mclock(r ; )) f 



so the average amount at the end of the day's business is just 



IP; 



/ L 



(iii) Average amount in store, (This is the 'true 1 average, avcragw Tn real time.) 
This is calculated simply as 

£"{a(clock(j)) + b{c\ock{s)))dsIL , 
again a candidate for the numerical integration routine. 
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Solution: Program Specifications 

There follows a largely complete Pascal program to carry out the calculations laid out in 
the mathematical description. 

1. Variables 

type 

vec «= array [0..MAXJ of real; 

daily * array [0..365] of real; 

vecjnt « array [0..M AX] of integer, 
(the parameter MAX should be set large enough to cope with the biggest value of nstock, 
the number of stockings. MAX * 700 will probably do. } 



(There follows a list of the variables used in the program, together with their type, range 
of possible values, and symbol used in the mathematical description (if applicable). It is 
assumed that all the variables above the line have already been input. } 



Variable name 


type 


range 


math, symbol 


nday 


integer 


1„„,365 


L 


openiimcs 


daily 


array of increasing reals in [0,1] 


iSj) 


closetimes 


daily 


array of increasing reals in [0,1] 


«y> 


hisunovement 


xeal 


positive integer 


n' 


CSL 


real 


(0,1) 


I* 


hisuCSI. 


real 


(0.1) 


a 


proportions 


daily 


reals in [0,1] 


<*•,) 


variability 


real 


positive real 


V 


nstock 


integer 


0„,.,MAX 


N 


stock.true.timc 0 ^ 


vec 


reals in [0,1] 


in) 


ndel 


integer 


0 MAX 


M 


del_tme_timc (1) 


vcc 


reals in [0,1] 


(D' t ) 


order.tme - rime <l) 


vec 


Ttalsin[0,l] 


TO 


delay_protect 


real 


[0,1] 


8 




real 


positive real 


C* 


target_stockJcvels (4) 


positive teals 


(wp 
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demand 


real 


positive real 




sigma 


real 


positive real 


0 


stock_shop_time 


vec 


reals in [0,1] 




dcI_shop_time 


vec 


reals in [0,1] 




order_shop_time 


vec 


reals in [0,1] 




D 
r 


vec 


reals in [0,1] 


(Pi) 


tnod_dcmands 


vec 






sales 


vec 


positive reals 


(Si) 


raw_tsl 


vec 


positive reals 




target_stockJevels 


vec 


positive reals 


(W.) 


effectives! 


vec 


positive reals 
• 


(Wj) 


shelf 


vec 


positive reals 


(Oi) 


backroom 


vec 


positive reals 




order.to 


vec 


positive reals 




var 


Teal 


positive real 


a 2 


h 


real 


positive real 




flag 


vecjnt 


-1,0,.. M MAX-1 




lcad_time_salcs 


vcc 


positive real 




lost_sales 


real 


positive real 




CSL_achieved 


real 


real in [0,1] 





0) The entries stock.true.timelO] stock_true_time [nstock-l] are the succes- 
sive times at which stockings are made. All subsequent entries in the array are Ini- 
tialised to 0. Similarly for dcLtruejime, order_true_time. 
w The program contains a facility for calculating the array stock jniejime from 
the input value stock_true_iime[0| alone; so in some uses, it would be satisfactory 
to input only this one value. 

There may be no value of C offered, because the main program has the capaci- 
ty to calculate C . , v ■ * i 
' 4) The user may wish to input the target stock levels (one for each of the ndel 
deliveries), but the program will calculate these if they are riot input 
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2. Initialisation 
begin 

epsilon :» 1.0E-6; 

variability := max (variability, epsilon); 
demand := hist_movcmeni/hist_CSL; 
sigma dcmand*variability; 
var sigma*sigraa; 
h :- 1.0/nstock; 

stock_true_tirac [nstock] > stock_Hue_time[0); 
del_twe_iime [ndel] := deLtrue_iime[0]; 
order_twe_dme (ndelj := order_true_time[Oj; 

(These conventions make 'wrap round' easier. } 

for j:»0 to ndel do 

begin 

del_shop_timcrj| :» clock(deUtrue_timcQ], opentimes, closetimes, 

proportions, nday); 
order«^hop_timeO] clock(oider_true_timelj]. opentimes, closetimes. 

proportions, nday); 

end; 

{ It may be that the user has input only stock_irue_timelOJ, and is leaving it to the program 
to calculate the rest of the array. This is done as follows: 

stock_shop_time[0] := clock(stock_troejime[0], opentimes, closetimes, 

proportions, nday); 

for j 1 to nstock do 
begin 

stock_shop_timcljl :« plus(stock_$hopjimeU-l].h); 
Stock_true.time[j] :» clockinv(stock_shop_timc[j1, opentimes, 

closetimes, proponions, nday); 1 

end; 

We assume therefore that the arrays stock_tructime and stock.shop.time have been ini- 
tialised.) 



for j :« 0 to nstock- 1 do 
begin 

plj] :- plus(stock_shop_tinjeO+l], -stock_.shop_time[fl); 
mod.demands(jl :■ modme*n(pfj1, tigma, demand) 
end; 
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Matn Routines 

3. Calculating shelf capacity 

{Firstly, note that if the shelf capacity is C, then the expected tost sales for the whole 
study period (assuming no backroom interference) is given by the function '* 
tot_lo$s(C.Teal;sigma:real; nstock:integer, p:vec; mod^demands:vec):reaL 

There are three methods available. 
Method 1 

Use a root-finding routine to pick C so that 
totJoss(C,...) - demand* (1-CSL). 

Method 2 

Find which of p[01,...,p[nstock-ll is greatest (p[rj, say) and then choose C so that 
loss(Csigma,p[rJ, mod_demands[r]) ■ demand*ptr]*(l-CSL). 

Method 3 

This is used when a retailer requires to be covered against the busiest consecutive K days 
trading 

read(K); 

x:«0.0; *•* 
fori:=0tonday-ldo 

begin 

y:-0.f>, 

fork:-0toK-ldo 

begin r> ' 

1 := (i+k) mod nday, 
y := y + proportions 0] 

end; 
if(y>x)ihen 

x:-y; 

end; 

z := modmean(x, sigma, demand); 
{and now we use the root-finding routine to pick Cso that 
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lo$$(C Sigma, x, z) = demand* x * (1-CSL). } 

(By now, the program has a value of C to work with, either calculated by one of the above 
methods, or input by the user. Now we can compute the expected sales in each interstock- 
ing period:) 

for j := 0 to nstock-1 do 

sale$[j] :■ dcmand*p(j] - lo$$(C Sigma, pfjK nK>djtemands[j]); 

4. Calculating target stock levels, and order-to levels 
for j := 0 to ndcl do 
begin 

deLtrue.timeO} plus(deLtrue_tkneIjJ, delay_protect -epsilon); 
deL$hop_time£j] ctock(del.trae_time[j], opentiraes, closetimes, 

proportions, nday) 

end; 

{next, we work out the effective delivery times, and store these in deLtnie_timc, 
del_$hop_time.} 

fori:=OtoMAXdo 

flagp] :* -1; 

fork:»Otondel-l do 
begin 

forj :«0tonstock-t do 
begin 

l:«(j+l)modnstock; 

if contain(deLtrue_time[k], stock.true.timelj], stock jrue-tijnjeOl) then 
begin 

deLtrue_time[k] :«plus(stock_true_timep] t -epsilon); 
deLshop_time[k] :*plus($tock_shopjimell] t -epsilon); 
flagdl :-k 

end; 

end; 

end; 

deljshop_time[ndel] :*dcL$hopjimel01; 
(now we get the order times into shop time) 



for k :« 0 10 ndel do . ■ ■ 

order_shop_time[k] :« pius(clock(oider_tme jim£[k], opentimes, cfosetimes, 

proportions, nday), -2* epsilon); 



[Now the calculation of target stock levels: } 

forj :»0tondcl*l do 
begin 

raw_tslIJl:»0.0; 

for i := 0 to nstock- 1 do 

begin 

if c*>ntain(stock_$hop_Ume[i] f dcl_shop_timeQ], deLshop^timefj+1]) 
then 

raw_tsl[j] raw_tsl[j] + salesli]; 

end; 

target_stockJevc!s[jJ :» maxOaw.tslljLC); 
(this assignment statement must be skipped if the user has chosen to input his own 
target_stockJevels} / 
end; 

(Now the calculation of lead-time sales: } 
fork:»OtondeMdo 
begin(l} 

x:=0.0; 

for j := 0 to nstock-1 do 
begin{2) 

if contain(stock_shopjime[j], order_shop_time[k] l deLshop_time[k]) 
then 

x:»x + sales[j|; 

else 

if . 

contain(order_shopjimelk],stock_shop_^ 

then 

bcgin{3) 
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t := plus(order_shop_time[k], •stock.shop.timelj]); 
z :■ loss(C*plj]/t, sigraa, pfj], mod_demands[j]); 
E:-i*(demand-fc/p[j])); 
z:=sales01 -z; 
x x+z; 



end;{3) 
end;{2} 

lead_time_sales[k] := x; 
end;(l} 

{ It may be that the user has input target stock levels, rather than relying on the values cal- 
culated by the program. Either way, there is an array target_stock_levels containing this 
data.) 

for k:=0 to ndel-1 do 

order_to[k] :» lead_time_sales[k] + targetjstock Jevcls[kl; 

5. Lost sales 

(Firstly, calculate effective.tsl.) 

x:=0.0; 
k:-0; 

for j :■ 0 to ndel-1 do 
begin 

if (targct_stock_levels[j] > x) then 



effectivc_tsl[kl := targetjuockjevelslk]; 

fori :■ 1 to ndel-1 do 

begin 

j := (k+i-1) mod ndel; 
1 := (k+i) mod ndel; 

x :«=effectivc_t$HjJ -min(raw_tslfjl, effective_tsHjl); 
effective.tsip] :- max(x, target_stockJevels[l]) 

end; 



begin 



k:-j; 

x :* target,stock_levels|j1 



end; 



end; 
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{From this, we calculate next the overall lost sales, and the customer service level 
achieved.) 

x:=0.0; 

fork:=Otondel-l do 

x min(raw„tsl[lc], effcctive.tsllkl) ♦ x; 
lostjsales :» demand • x; 
CSL.achieved := x/demand; 

6. Average inventory calculations 

fork:»Otondel-l do 
begin 

for j :» 0 to nstock- 1 do 

begin . 

if(flagQ]«k)then 
begin 

backroomOJ := rnax(effective_tslIki-C, 0.0); 
shelflj] min(C,effective_tsipc]) 

end; 

end; 

end; 

{this has set the values of backroom and shelf at stocking times when there is a delivery. 
The next task is to calculate U for all stocking times, for which wc must find the first 
stocking rime which has a delivery, arid work from there.) 

i:=0; 

while (flag[il--l) do 

i:«i+l; r> • 

for j := 1 to nstock* I do 
begin 

k := 0+j- 1 ) mod nstock; 
1 := (i+j) mod nstock; 
if(flag[l]«-l)then 
begin 

x min(shelffkl, salesM); 
y := shelflkj - x; 

backroomni rnax(backxoom[kj - x, ttO); 
shelflU := min(backroom[k] + y, C); 
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end; 

. end; 

(The functions a and b allow us to obtain for any t in [0.1] the average amount on the 
shelf (respectively, in the backroom) at time t 

To calculate the average amount seen by the customers, we do a (wpcziuowule) calcula- 
tion of Pa (0<a. 

To calculate the average amount in the store at the end of the day's business, we take 

forj:«0tonday-ldo 
gfj] a(clock(j/nday)) + b(clock(j/nday)); 
average := 0.0; 
forj :»0tonday-ido 
average := average + gO I; 
average := average/nday, 
To find the genuine time-average, we make the integral 
^(a(clock(t))+b(clock(t)))dt } 
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Subprogram Specifications 

function clock(t: real; a: daily; b: daily; pr. daily; n; integer): real; 

{given a time t in the interval (0,1), this function calculates the proportion taf business 

done by time t. It also extends periodically to the whole line.) v 

var 

mo : integer, 
x t sum, y : real; 
begin 

m:=trunc(t); 
x:»t-m; 

if(x>=b[n-l]) then 
clock:- 1.0 

else 
begin 



y :» max(0.0, (x - a01V(blj] - »0D); 
clock := sum + prfjl*min(1.0,y) 

end; 

clock := clock + m; 

end; 

We also need 

function max(x; real; y: real): real 
which returns the larger of x and y, 

function min(x: real; y: real): real 
which returns the smaller of x and y, and 



j:-0; 

sum :- 0.0; 

while (x>=b01) do 



begin 



sunt:* sum + pr£jl; 
j J+l 



end; 
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function frac(x: real): real 
which returns x * trunc(x), the fractional part of x. 



function plus(x: real; y: real): real; 
(this adds x toy modi) 
begin 

plus frac(x+y) 

end; 



function clockinv(t: real; a: daily; b: daily, pr: daily; n: integer): real; 

{this is the inverse function to clock; given t, it returns $ such that clock(s) • t The array 
a contains opening times, b contains closing times, and pr contains proportions of busi- 
ness done on the various days) 

var 

mj: integer, 
x, sum: real; 
begin 

m :» trunc(t); 
x :» t-m; 



end; 
z:» sum -x; 

clockinv :» m + b[j] - (b[j] * a[j])*z/pr[}]; 

end; 

function N(z: real): real; 

(this gives a good approximation to the integral from z to infinity of the standard nonnal 
distribution.) 



j:«0; 

sum :» 0.0; 

while (x >-sum)do 



begin 

sum :» sura+pr[j]; 
j:=j+l 
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var . . : 

x t y,v:real; 
begin 

y:-l/(l +0.23164190*abs(z)); 
x :- 0.3989423*exp(-0.5»z*z); 

v:-u30274*y+y*y*y*y; 

v:«v-1.821256*yVy*yi 
v:»v+ 1.78U78*y*y*yi 
v:-v-0.356538*y*y. 
v:-v + 0.3!93875*y; 

v:»x*v; * 
if (z>* 0.0) then 

N:»v 

else 

N:-l-v 

end; 

function A(z: real): real 

{this gives the expectation of the positive part of Y-t, where Y is a standard normal ran- 
dom variable) 
begin 

A :=» 0.3989423*exp(-0.5*z*z) - z*N(z) ; : '' 

end; 

function modmean(x: Teal; s: real; d: real): real; \ 
(if x is the proportion of the business done in some imerstocking interval, s is the, standard 
deviation of the overall demand, and d is the overall mean demand, this function works 
out the mean of a normal random variable Y so tba» the expected value of Y* is equal to 
xd.l 

Once again a general root-finding touting is needed. This function should return that 
value y such that 

• 

Hie value y will always be less than xd; negative values are possible, 
function loss(C: real; sigma: real; x: real; mean: real): real; 
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{this calculates expected lost sales in our intcrstocking interval when the shelf capacity is 
C and a proportion x of the business for whole study period is done} 

var 

z: real; 
begin 

2 :« sigma*sqrt(x); 
loss :« 2* A((Gtnean)/z) 

end; 

function totJoss(C: teal; sigroa: real; nstock: integer, p:vec; mod_mcan: vec): real; 
var 

j: integer, 
begin 

toiJosssO.O; 
forj:«0tonstock-ldo 

toUoss :» toUoss + lossCC, sigma, p[fl, mcxLmean[jl) 

end; 

function contain(x: real; a: real; b: teal): boolean; 

(this function should only be used on arguments in the interval [0,1]; it returns 'true' if x 
is in the interval from a to b (with the endpoints of the interval wrapped around) other- 
wise 4 f*lse\} 
begin 

if(a<b)then 

if(a<-x)and(x<*b)then . 

contain := true 
else contain false; 

else 

if (a<=x)or(x<=b)thcn 

contain :« true 
else contain false 

end; 

function aft: real; stock^shopjime: vec; shelf: vec; sales: vec; demand: real): real; 
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{ this works out the average amount on the shelf at time t ) 
var 

j: integer, 
begin 

forj:«Otonstock-l do 
begin 

If contain^ stock_shop_time[j] f stock_shopjimc|j+l]) then 



a := roax($helf[j] - demand*(t-stock_shop_timeOt), max($helflj]*$alcs[jl, 



function b(t: real; stock_shop_time: vec; backroom: vec): real; 
(this works out the average amount in the backroom at time t) 
var 

j: integer, 
begin 

for j :« 0 to nstocH do 
begin 



if contain^ $toc!cjhopjime[j}, stod^shop_time|j+l]) then 
b:»backioom(j] 



0.0)) 



end; 



end; 



end; 



end; 
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While the invention has been described with ref- 
erence to details of the illustrated embodiment, these de- 
tails are not intended to limit the scope of the invention 
as defined in the appended claims. 



? 



4 
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Claims 

1 1. A retail shelf inventory management system 

2 comprising: 

3 memory means for storing inventory .data, ,said in- 

4 ventory data including a study period and a retail opening 

5 time and a retail closing time for each day within said 

6 study period and historical inventory data, for sai<3 study 

7 period for each of a plurality of inventory items including 

8 a movement value, a customer service level, a proportion of 

9 business done for each day within said study period, and a 

10 variability of demand; and 

11 optimizing means operatively connected; to said 

12 memory means for optimizing a shelf capacity for each of 

13 said plurality of inventory items. 

1 2. A retail shelf inventory managem^git system 

2 as recited in claim 1 wherein said optimizing meai^ns include 

3 means for calculating an expected lost sales value for said 
J 4 study period and identifying said shelf capacity resppnsive 

5 to said calculated expected lost sales value. 

. 

1 3. A retail shelf inventory management system 

2 as recited in claim 1 wherein said inventory data includes 

3 a desired customer service level value to be achieved; and 

4 wherein said optimizing means include means for calculating 

5 an expected lost sales value for said study period and 

6 identifying said shelf capacity responsive to said calcu- 

7 lated expected lost sales value and said desired customer 

8 service level value. 
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1 4. A retail shelf inventory management system 

2 as recited in claim 1 wherein said inventory data includes 

3 stocking time values and wherein said optimizing means 

4 include means for identifying an interstocking period 

5 having a maximum proportional value of business done and 

6 . identifying said shelf capacity responsive to said 

7 identified interstocking period. 

1 5. A retail shelf inventory management system 

2 as recited in claim 1 wherein said inventory data includes 

3 a desired customer service level value to be achieved and 

4 said historical inventory data includes interstocking time 

5 values and wherein said optimizing means include means for 

6 identifying an interstocking period having a maximum pro- 

7 portional value of business done and identifying said shelf 

8 capacity responsive to said identified interstocking period 

9 and said desired customer service level value. 

1 6. A retail shelf inventory management system 

2 as recited in claim 1 wherein said inventory data includes 

3 a desired customer service level value to be achieved; and 

4 wherein said optimizing means include means for identifying 

5 a predetermined number of consecutive days having a maximum 

6 proportional value of business done and identifying said 

7 shelf capacity responsive to said identified consecutive 

8 days and said desired customer service level value. 
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1 7. A retail shelf inventory management system 

2 as recited in claim 1 wherein said inventory data includes 

3 a number of shelf stockings and at least a first stocking 

4 time; and further comprising means for qalculating' an' 

5 optimum shelf stocking schedule for said study ""period. 

1 8. A retail shelf inventory management system 

2 as recited in claim 1 wherein said historical inventory 

3 data includes a number of order deliveries and delivery 

4 times in said study period for each of said plurality of 

5 inventory items and further comprising means for caiculat- 

6 ing an effective delivery time and means for calculating a 

7 target stock level responsive to said calculated effective 

8 delivery time and said shelf capacity fbr any of said 

9 inventory items. 

1 9. A retail shelf inventory management system 

2 as recited in claim 1 further comprising means for calcu- 

3 lating an order to level defining an order quantity amount 

4 for any of said inventory items. 

1 10. A retail shelf inventory management system 

2 as recited in claim 1 further comprising means for calcu- 

3 lating an expected lost sales value responsive to a user 

4 selected shelf capacity. 

1 11. A retail shelf inventory management system 

2 as recited in claim 1 further comprising means for calcu- 

3 lating an expected lost sales value responsive to said op- 

4 timized shelf capacity. 
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1 12. A retail shelf inventory management system 

2 as recited in claim 1 further comprising means for calcu- 

3 lating an average shelf inventory value responsive to said 

4 optimized shelf capacity. 

1 13. A retail shelf inventory management system 

2 as recited in claim 1 further comprising means for calcu- 

3 lating an average store inventory value responsive to said 

4 optimized shelf capacity, 

1 14. A retail shelf inventory management system 

2 as recited in claim l further comprising means for calcu- 

3 lating an average shelf inventory value and an average 

4 store inventory value responsive to a user selected shelf 

5 capacity. 

1 15. A retail shelf inventory management system 

2 as recited in claim 1 further comprising input means cou- 

3 pled to said memory means for receiving user input selec- 

4 tions. 

1 16. A retail shelf inventory management system 

2 as recited in claim 1 further comprising display means cou- 

3 pled to said optimizing means for displaying results for a 

4 user of the system. 
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1 17. A retail shelf inventory management system 

2 comprising: 

3 memory means for storing inventory data, said in- 

4 ventory data including a study period .and a retail opening 

5 time and a retail closing time for each day within said 

6 . study period and historical inventory data for said study 

7 period for each of a plurality of inventory items including 

8 a movement value, a customer service level, a prepdrtion of 

9 business done for each day within said study period, and a 
10 variability of demand, a number of deliveries and order de- 
ll livery times within said study period and a number of 

12 stockings and stocking times within said study period, and 

13 a customer service level to be achieved; 

14 optimizing means operatively connected to said 

15 memory means for optimizing a shelf capacity for any of 

16 said plurality of inventory items; and 

17 means coupled to said optimizing means for calcu- 

18 lating an effective delivery time and an order quantity 

19 amount for any of said inventory items responsive to said 

20 optimized shelf capacity. 

1 18. A retail shelf inventory management system 

2 as recited in claim 17 wherein said optimizing means in- 

3 elude means for calculating an expected lost sales value 

4 for said study period and identifying said shelf capacity 

5 responsive to said calculated expected lost sales value. 

1 19. A retail shelf inventory management system 

2 as recited in claim 17 wherein said inventory data includes 

3 a desired customer service level value to be achieyed; .and 

4 wherein said optimizing means include means for calculating 

5 an expected lost sales value for said study period and 

6 identifying said shelf capacity responsive to said caicu- 

7 lated expected lost sales value and said desired : customer 

8 service level value. 
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1 20. A retail shelf inventory management system 

2 as recited in claim 17 wherein said inventory data includes 

3 a desired customer service level value to be achieved and 

4 said historical inventory data includes interstocking time 

5 values and wherein said optimizing means include means for 

6 • identifying an interstocking period having a maximum pro- 

7 portional value of business done and identifying said shelf 

8 capacity responsive to said identified interstocking period 

9 and said desired customer service level value. 
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